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ABSTRACT 

We discuss the determination of the Unear power spectrum of dark matter fluctuations from measure- 
ments of the transmitted flux Ly-a Forest power spectrum. We show that on most scales probed by 
current measurements, the flux power spectrum is significantly affected by non-linear corrections to the 
linear dark matter power spectrum due to gravitational clustering. Inferring the linear dark matter power 
spectrum shape is therefore difficult due to non-linear corrections driving the non-linear power spectrum 
to a k~^''^ shape nearly independent of initial conditions. We argue that some methods used in previous 
estimates underestimate the uncertainties in the shape of the linear dark matter power spectrum. 

Subject headings: cosmology: theory - intergalactic medium - large scale structure of universe; quasars 
- absorption lines 



1. INTRODUCTION 

In the last few years the Lyman a forest has become 
increasingly important as a probe for cosmology. The 
availability of high resolution high signal-to-noise spectra, 
mainly from the HIRES spectrograph on the Keck tele- 
scope, has revolutionized the field. The data have been 
compared with the theoretical predictions of currently fa- 
vored cosmological models, mainly using numerical sim- 
ulations (Cen et al. 1994, Hernquist et al. 1995, Zhang 
et al. 1995, Miralda-Escude et al. 1996, Muecket et al. 
1996, Wadsley & Bond 1996, Theuns et al. 1998). The 
success of these models has been impressive. The ob- 
served properties of the forest are very well reproduced 
even though the models were primarily constructed in the 
context of galaxy clustering - this was a non-trivial accom- 
plishment. Examples of comparisons of model predictions 
with observational data can be found in Croft et al. (1998), 
Bryan et al. (1999), McDonald et al. (1999, 2000), Nusser 
& Haehnelt (2000), Schaye et al. (1999, 2000), Croft et 
al. 2000, Meiksin, Bryan & Machacek (2001), Choud- 
hury, Srianand & Padmanabhan (2001), Zaldarriaga, Hui 
& Tegmark (2000) (ZHT2000), and Pichon et al. (2001). 
The same cosmological models also produce acceptable fits 
to many other cosmological observations such as the power 
spectra of CMB fluctuations, the luminosity distance to 
high redshift supernovae, the number density of massive 
clusters, etc. (see Wang, Tegmark & Zaldarriaga 2000 for 
a recent summary). 

The success of the numerical models of the forest has 
also led to the construction of several analytic approxima- 
tions (e.g. Bi & Davidson 1997, Gnedin & Hui 1996, Hui 
et al. 1997b, Viel et al. 2001, Schaye 2001). The physics 
that is responsible for the absorption seems to be well un- 
derstood so the hope has emerged that we could use the 
Lyman-a forest to further constrain cosmological models. 
The advantages of the forest over other observations are 
clear, since we get to observe the universe at an earlier time 
we may be able to infer the rate of growth of density fluctu- 
ations. It is also assumed that because we are observing at 
a higher redshift, scales that today (e.g. in galaxy surveys) 



are in the "non-linear regime" , can be probed when they 
are still "linear", thus greatly simplifying the analysis. In 
this paper, we challenge this assumption. 

Studies of the forest have other applications. Through 
the Alcock-Paczynski test (Alcock & Paczynski 1979) ap- 
plied to the correlation functions for pairs of lines of sight, 
it seems possible to directly infer the geometry of the uni- 
verse (Hui, Stebbins & Buries 2000, McDonald & Miralda- 
Escude 2000, McDonald 2001). Higher order statistics of 
the flux may also prove useful in testing that the structure 
observed in the flux is actually produced by gravitational 
instability (Zaldarriaga, Seljak & Hui 2001) and in testing 
if there are fluctuations in the temperature of the IGM, 
left for example as remnants of a late ionization of Helium 
II (Zaldarriaga 2001, Theuns et al. 2001). 

The power spectrum of the flux is perhaps the most 
studied of the statistics of the Lyman a forest flux. In- 
ferring cosmological parameters from the power spectrum 
is however a difficult challenge. A naive approach would 
be to run hydrodynamical simulation for enough models 
to cover the parameter space of interest and compare the 
results with the observed data. This approach is imprac- 
tical as the computer time needed to run the number of 
simulations required is enormous. One should also remem- 
ber that on top of the usual cosmological parameters that 
are needed to specify the cosmological model {i.e. the en- 
ergy densities in the different components and the power 
spectrum of initial density fluctuations) in the case of the 
forest one needs to specify the time evolution of the back- 
ground of ionizing radiation that is responsible for ionizing 
and heating the gas. This requires the introduction of ad- 
ditional parameters. One should also keep in mind that 
even if we were able to run a large enough grid of hydrody- 
namical simulations, it is not obvious that all the relevant 
physics is correctly included in these simulations, "feed- 
back" effects from the forming stars and galaxies being 
the obvious example (see e.g. Meiksin, Bryan & Machacek 
2001 for possible conflicts with observations). 

A significant simplification can be achieved with the in- 
troduction of the N-body model of the forest. In this 
model the spatial distribution of the dark matter is ob- 
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taincd by running an N-body simulation. The gas is then 
assumed to traee the dark matter (or maybe a smootlied 
version of the dark matter field to account for the effect 
of gas pressure). Furthermore it is assumed that the gas 
is in photo-ionization equilibrium and follows a power-law 
equation of state. This simple prescriptions lead to a sim- 
ulated Lyman a forest that is in reasonable agreement 
with the predictions of full hydro simulations. As could 
have been expected, the agreement is good but not per- 
fect (Meiksin & White 2000). The great advantage of the 
N-body approximation is that it significantly reduces the 
computational burden of exploring parameter space. A 
first attempt at implementing the straightforward method 
of running a grid of N-body simulations and comparing 
that with the forest data to extract parameters was pre- 
sented in ZHT2000. 

In the quest of incorporating data from the forest to 
studies of cosmological parameters even the N-body model 
is computationally too expensive. This has lead to the 
introduction of approximate techniques to invert the ob- 
served flux power spectrum to get the 3D linear matter 
power spectrum of the dark matter. The results of these 
inversions are then compared to any model of interest and 
the results of these comparisons are used to constrain pa- 
rameters. It is clear that the inversion method greatly sim- 
plifies the comparison between theory and observations, 
which only requires the ability to predict the linear power 
spectrum of any given model, a straightforward task. Ob- 
viously the relation between the linear dark matter power 
spectrum and the flux power spectrum is non-trivial, so 
these inversions are not rigorous inversions in the math- 
ematical sense. The resulting prescriptions, however, are 
calibrated with numerical simulations, and are the result 
of detailed studies of the physics of the forest (Croft et al. 
2000). 

The aim of this paper is to point out the limitations of 
inversion techniques which attempt to match predictions 
of the forest from knowledge of the linear power spectrum 
of dark matter. In particular, we point out that the evo- 
lution of fluctuations at 2; ~ 3 is significantly non-linear 
at scales that can be orders of magnitude larger than esti- 
mated by naive linear arguments. Furthermore, non-linear 
corrections tend to drive the density (and velocity) field 
to "universal" spectra nearly independent of initial condi- 
tions, erasing information and leading to weak constraints 
on the shape of the primordial dark matter power spec- 
trum. 

The rest of this paper is organized as follows. In §2 we 
summarize the constraints on the primordial power spec- 
trum that we obtain by running a grid of N-body simula- 
tions and compare them with those obtained by the inver- 
sion technique described above. The remaining sections 
present what we believe to be the main reasons behind 
these results. In §3 we discuss the role of non-linear cor- 
rections, in §4 we discuss the effects of the transformation 
between density and flux and the importance of redshift 
distortions. We conclude in §5. 

2. LIKELIHOOD FROM N-BODY SIMULATIONS 

Our first objective is to convince the reader that there 
are problems with the error bars that result from inversion 
methods. We use the technique described in ZHT2000 to 
find constraints on the amplitude and slope of the linear 



power spectrum of dark matter fluctuations. Basically we 
run a grid of P'^M N-body simulations of the standard 
CDM model in boxes of size L = 16 /i~^Mpc. We choose 
a range of spectral indices for the primordial power spec- 
trum, (rip = 0.4,0.7,0.9,1.0,1.1,1.3,1.7), rip = 1 corre- 
sponding to the usual scale-invariant spectrum, and store 
the outputs at several different times, corresponding to 
scale factors a = (.05, .065, .085, .11, .14, .19, .24, .31, .41). 

The relation between gas density and dark matter den- 
sity is modeled by introducing additional parameters. 
To model the effects of pressure we smooth the density 
field with a three-dimensional smoothing scale we call 
kf. For simplicity we used a Gaussian window, W{k) = 
exp{-{k/kff as described in ZHT2000. The smoothing 
scale kf is related to the Jeans scale kj by, 



kf = rjkj; kj = ac^ ^ \J ^ixGp, 



(1) 



where Cg is the sound speed of the gas and 77 is a function of 
time and the details of the ionization history with values 
around 77 ~ 2 for standard assumptions (Gnedin & Hui 
1998). Equation (1) leads to. 
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(2) 



where rup is the mass of the proton. 

The equation of state of the gas is described by two pa- 
rameters, the temperature at zero overdensity Tq and the 
slope a, T = ro(l -I- (5)" where i5 is the overdensity (Hui & 
Gnedin 1997a). The optical depth of a given fluid element 
is given by r = ^(1 -\- 6)^, the assumption of ionization 
equilibrium leads to /3 = 2 — 0.7a. The amplitude A, which 
physically depends on the strength of the ionizing back- 
ground, is chosen so that the constructed spectra have the 
observed mean transmission which we call F. For each 
simulation we construct grids of model predictions with 
{kf ,Tq, a, F) as parameters. Our full parameter vector 
p = (a. Up, kf,To, a, F) has 6 dimensions. 

For each model in our six dimensional space we calcu- 
late the predicted flux power spectrum P/(p). We com- 
pare these prediction to the data and compute a likelihood 
function as described in ZHT2000. We used a Gaussian 
approximation for the likelihood, 



£(d; p) oc exp 



di - Pfiip) 



(3) 



where i runs over the different data points di, the model 
prediction for that wavevector are -P/i(p) and ai are the 
error bars on each point. To obtain one dimensional con- 
straints on individual parameters or two dimensional con- 
straints for a pair of parameters, we marginalized the like- 
lihood one parameter at a time, reducing the dimensions 
of the grid by one in each step until only the parameter or 
pair of parameters of interest was left. 

When marginalizing over parameters we allow F to 
vary within observational constraints. At z = 2.72 cor- 
responding to the mean redshift of the sample used we 
take F = 0.707 ± 0.035, a 5% uncertainty. We should 
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note that wo have added in quadrature a "theoretical un- 
certainty" to the observational errors. This extra uncer- 
tainty was estimated from the fluctuations in the mean 
flux computed in different realizations of the same model 
and could be improved upon with better simulations. We 
include the constraints on F by multiplying the likelihood 
coming from tjie forest power spectrum by exp(— 5x^/2), 
with Sx^ = [{F - 0.707)/0.035]2. 

In ZHT2000 we presented constraints on the primordial 
power spectrum. In that paper we were mainly interested 
in the effects of the temperature and the constraints that 
could be obtained on it. The temperature smooths the 
spectra and we wanted to make sure that we did not un- 
derestimate the uncertainties by assuming we knew what 
kf, the other source of smoothing, actually was. For that 
reason we had a grid in kf that was too conservative, 
kf = (5,10,20,30,40,45,50,55,60,70,80)/i Mpc-^ The 
smallest values in this grid (the largest smoothings) pro- 
duce much more smoothing than can be expected from 
pressure. Howevcir because kj smooths the field in 3D 
before we apply the non-linear transformation that takes 
density to optical depth and flux, while Tq acts in ID 
smoothing the optical depth, we were able to distinguish 
between them. We also found that the kf smoothing was 
quite degenerate with the spectral index of the perturba- 
tions n so our constraints on that parameters were weak. 

In this paper we improve on the earlier analysis by using 
newer data from Croft et al. (2000). We use the power 
spectra of their "fiducial sample" which has mean red- 
shift z = 2.72 (table 3 of Croft et al. 2000). The sample 
was compiled using data from the Low Resolution Imaging 
Spectrometer (LRIS) and the High Resolution Spectrom- 
eter (HIRES) on the Keck Telescope. Data in the redshift 
range 2.3 < 2; < 2.9 with a total path length 5z « 25 was 
used to create the "fiducial sample" (Croft et al. 2000). As 
another improvement in our marginalizations we restrict 
the kf grid to values kf > 30h Mpc^^. With the improved 
data and a more physical "prior" for kf we hope to find 
tighter constraints on the primordial power spectrum. 

In figure 1 we show our constraints in the — Ucs, the 
amplitude and slope of the linear power at a scale kp = 
0.03 (km/s)-i. That is A'^{k) = Ank^Pik) and /^^{kp) = 
Al with ((iln A^(fc)/rflnfc)|fcp = ne// + 3. To illustrate the 
effect of the uncertainty in A:/ we show the results if we 
marginalize over kf with the restriction fc/ > 30h Mpc~ 
or if we just assume kf = 30/i Mpc~^. For comparison we 
also plot the results obtained by Croft et al. 2000 with 
the inversion technique. Note that when we marginalize 
over kf, that is when we allow smaller kf smoothing scales, 
models with smaller spectral indices become allowed. This 
makes sense because one can compensate for the lack of 
power in the dark matter by assuming that the smoothing 
due to pressure is smaller. 

Figure 1 shows that we find acceptable models outside 
the range inferred from the inversion technique. In figure 
2 we show examples of such models. The two models have 
very different spectral indices of the linear power spec- 
trum, one has Up — 0.7 (rioff = —2.7, = 0.42) and the 
other has Up = 1.3 {ues = — 2.3, A^ = 0.57). Part of the 
increase in the allowed range of A J seen in figure 1 comes 
from our more conservative choice of errors for F. Unfor- 
tunately, the relation between the amplitudes of the dark 



matter and flux power spectrum is very sensitive to F, a 
quantity that is sensitive to small scale physics and may 
not be very well modeled by dark matter only simulations. 
Figures 1 and 2 clearly show that the inversion technique 
underestimates the errors, thus results from such studies 
should be interpreted with care if one wants to incorporate 
them to likelihood analysis of cosmological parameters. 

It is important to understand that the non-linear cor- 
rections make the power spectra of the two models under 
consideration look similar but that in order to fit the data 
we must allow some extra freedom in another parameter 
to accommodate the residual differences. In the example 
shown, a had to vary between 0.1 and 0.3. This means 
that if there are independent ways (not using the power 
spectrum) to measure a, and in general the equation of 
state, then the error-bar on n can conceivably be reduced. 
Using the b-parameter-column density distribution might 
be one way of achieving this, but the error on a is at 
present still quite large (eg. Ricotti et al. 2000, Schaye et 
al. 2000, McDonald et al. 1999, Bryan et al. 1999, Hui & 
Rutledge 1999). 

We should emphasize that the error-bars we find using 
our grid method are likely to underestimate the real un- 
certainties because we do not include uncertainties associ- 
ated with hydrodynamic effects. Nonetheless, it is possible 
that our uncertainties are larger than they should be. For 
example, the smoothing due to pressure (fc/), the temper- 
ature of the gas and the slope of the equation of state 
can all be determined once the ionization history is known 
[i.e. equations (1) and (2) for kf]. However we are allowing 
them to vary independently in our marginalization. This 
means that if we were to factor in some "physical pri- 
ors" that restricted these parameters or forced relations 
between them, our constraints could get tighter. 

The aim of the following sections is to try to identify 
what is causing the increased uncertainties. 

3. NON-LINEAR CORRECTIONS: CRITICAL INDEX 

The scales probed by the Lyman-a forest at 2; ~ 3 corre- 
spond to values of the linear mass power spectrum Phik) 
where 47rfc'^Pi(fc) ^ 0.01 — 1. This may suggest that these 
scales are safely in the linear regime. However, the am- 
plitude of linear fluctuations is not the only quantity that 
determines the magnitude of non-linear corrections. The 
other important quantity is the spectral index at scales 
in the weakly non-linear regime (Makino, Sasaki & Suto 
1992; Lokas et al. 1996; Scoccimarro & Frieman 1996). 
For typical CDM models, the effective spectral index rieff 
at the non-linear scale (fcni such that 47rfcJ5jPL(fcni) = 1) 
is quite negative riofr ^ —2.5 (e.g. compared to z = 
where n^s ^ —1.5). This means that non-linear correc- 
tions are large, since for steep spectra gravitational col- 
lapse is driven by large-scale coherent flows which substan- 
tially enhance the growth of density perturbations com- 
pared to the linear case. In fact, for power-law initial 
spectra we can write the non-linear mass power spectrum 
in one-loop perturbation theory (PT) as (Scoccimarro & 
Frieman 1996) 



P(fc) = P4fc)[l + a(n)(A) 



(4) 



where a{n) is a function of spectral index which mono- 




Fig. 1. — Constraints in the A — rieff plane. The solid line shows 
constraints when we &x. kf = 30h Mpc~ , the dashed curve is when 
we allow kf > 30ft Mpc~^. In both cases the lines correspond to 
an increase of Sx^ = 6.17 with respect to the best model (95% con- 
fidence for a Gaussian). The dashed circle are constraints obtained 
using the inversion technique and the same data by Croft et al. 
2000. 
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Fig. 2. — Transmission power spectrum from Croft et al. (2001) 
together with two models that are acceptable fits. The two models 
have parameter vectors p = (0.19,0.7,30,250,0.1,0.684) and p = 
(0.11, 1.1, 30, 250, 0.3, 0.684). 



power spectrum retains its linear shape. For n > —1.4, 
the non-hnear corrections are actually negative, so the 
non-linear power spectrum grows slower than in linear 
PT, becoming steeper than the linear spectrum. As a re- 
sult of this, the non- linear power spectrum is driven to- 
wards a "critical index" ric = —1.4, regardless of its linear 
slope (Scoccimarro & Frieman 1996). Although Eq. (4) 
is valid rigorously only for a power-law spectrum, CDM 
spectra behave in the same way, as seen in Fig. 3. Note 
how the shape of the non-linear power spectrum becomes 
~ fc~^'^ (dashed line) even at k < kni- In fact all the 
scales shown in the figure arc below k^i for the model with 
Hp = 0.7. Similar behavior is seen in Fig. 1 of Narayanan 
et al. (2000) and Fig. 3 of White k Croft (2000) for CDM 
models with truncated linear spectra. 

Figure 4 illustrates the importance of non-linear correc- 
tions as a function of power spectrum normalization cjg, 
for linear power spectra with shape parameter F = 0.21 
(dashed), F = 0.15 (dot-dashed) and F = 0.1 (solid), 
which approximately span the range of spectral indices 
obtained in figure 1. The calculation was done using one- 
loop perturbation theory^. The scales where linear (top 
three lines) and non-linear (middle lines) fluctuations are 
of order unity and where non-linear corrections are 20% 
(bottom lines) all agree at large as, where the spectral in- 
dex at the non- linear scale is about Ues ~ —1.5. At low erg, 
the effective spectral index at weakly non-linear scales be- 
comes much more negative, leading to stronger non-linear 
corrections. The wavenumber at which non-linear correc- 
tions are 20% can thus be a very small fraction of the 
non-linear scale (estimated either from the linear or the 
non-linear power spectrum). For the normalizations cor- 
responding to z ^ 3 (denoted by the horizontal bar in 
figure 4), the wavenumber at which non- linear corrections 
become important can be orders of magnitude smaller than 
that estimated from linear perturbation theory. In partic- 
ular, note that models with more negative spectral index 
(lower F) require higher erg to match the fluc;tuation am- 
plitude at the scales probed by the Ly-a forest, thus these 
models are most strongly affected by non-linear effects, at 
wavenumbers as small as A: ^ 0.004 km/sec. 

The outcome of this situation is that at z ^ 3 non-linear 
corrections are important at scales of interest for the Ly- 
a forest and drive the mass power spectrum to a k~^ '^ 
dependence. Therefore, measuring the Ly-a flux power 
spectrum at these scales docs not give a direct probe of 
the initial, linear, shape of the dark matter power spec- 
trum, leading to degeneracies in the range of initial shapes 
constrained by the data, as shown in figure 1. This is illus- 
trated rather clearly by comparing figures 2 and 3, where 
significantly different linear spectra display very similar 
transmission power spectra. However, before we can es- 
tablish that non-linear dynamics is the main reason behind 
the observed behavior, we must consider other non-linear 
effects that potentially enter: the non-linear mapping be- 
tween the dark matter density field and Ly-a transmission, 
and the effects of reclshift distortions. 



tonically decreases with n, becoming larger than unity for 
n < —1.7. For n < —1.4, the non-linear correction is pos- 
itive, leading to a non-linear spectrum less steep than the 
linear one. For n w —1.4, a becomes zero and thus the 



4. FLUX CORRELATIONS; A SIMPLE MODEL 

The discussion in the previous section was in terms of 
the dark matter density field 6. However, as we show be- 

^If anything, this calculation underestimates non-linear effects. 
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Fig. 3. — Linear (dotted) and non-linear (solid) dark matter 
power spectra for two models with primordial spectral index Up = 
0.7 (top) and rip = 1.3 bottom, corresponding to the transmission 
power spectra shown in figure 2. Note how non-linear corrections 
are very large for the Up = 0.7 case, even though the linear spectrum 
suggests linear PT should be accurate. Non-linear dynamics tends 
to establish a power spectrum with a critical index n = —1.4, as 
shown by the dashed line. 



low, there is a relevant range of scales for which the pic- 
ture described above applies as well to the transmitted 
Ly-a flux power spectrum. The reason is that the char- 
acteristic wavenumber for which the non-linear nature of 
the mapping between S and F enters, the smoothing scale 
wavenumber^ fcsM ^ 0.05 (km/scc)"-'^, is much larger than 
the wavenumber where non-linear corrections in the den- 
sity field become important (see figures 3-4). 

We take a simple toy model to illustrate our arguments. 
The transmitted flux F = exp(— t) is approximated as 
F = cxp[— A(l + (5)^] and S is the density field assumed to 
be Gaussian for the moment (we will relax this assumption 
below). In this case, the mean transmission is given by 



F = 



exp 



2A 



1 + 2Aa^ 



Wc work in terms of the flux contrast, Sp 
then the two-point function of dp reads, 



(5) 

F/F - 1, 



(1 + 2Aa^) 



(6) 



^1 -F4y4a2-h 4^2(^4 _ ^2) 

■^The sharp, exponential decline in the flux power spectrum seen 
in figure 2 is due to thermal broadening. The scale fcsM is related 
to the thermal broadening scale kj^jj ai 0.1{km/ sec)~^ by the non- 
linear mapping from optical depth to transmission power spectrum; 
this follows from Eq. (7) below. That is, for a given value of F and 
<T^, ^p/f^F is a function of i;/<J. Then, smoothing in the optical 
depth (approximately density) can be translated into smoothing in 
the transmission power spectrum. For cr > 1, ksM < kTH- 
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Fig. 4.- The importance of non-linear effects, as a function 
of (jg for linear power spectra with shape parameter F = 0.21 
(dashed), F = 0.15 (dot-dashed) and F = 0.1 (solid). The three 
lines denote the non-linear scale according to linear spectrum (top), 
non-linear spectrum (middle), and the scale at which non-linear 
corrections are 20%. At low redshift, high erg, all these scales co- 
incide, approximately. At high redshift, small erg, non- linear cor- 
rections arc much larger and thus become important at wavcnum- 
bcrs that can be a very small fraction of the non-linear scale. The 
horizontal bar denotes approximately the range of normalizations 
for 2 3. Note that translation from h/Mpc to inverse velocity 
units is given by k^ = k/(aH) Ri A:///o[f^m(l + z)]^^^'^ at high z, 
where = nm(z = 0). For ACDM with Q.'^ = 0.3, = 0.7, 
kv = fe/111.5 at z = 3. 



X exp 



4^2^ 



il + 2Aa^)[l + 2A{a^+^)] 



where ^ is the two-point function of the density field. Note 
that in the limit where ^ ^ (but not necessarily ^ < 1, 

since Acr"^ cannot be large due to the mean flux in Eq. (5) 
being of order unity), the above expression reduces to. 



4^2 



(1 + 2A<t2)2 



(7) 



Thus, in this case the flux correlation function is propor- 
tional to the density correlation function with some bias 
factor. This is certainly expected at large scales, since the 
flux is a local transformation of the density contrast, it fol- 
lows from well-known arguments (e.g. Fry & Gaztaiiaga 
1993, Scherrer & Weinberg 1998) that the correlation func- 
tion of the flux should be proportional to that of the un- 
derlying density field. The bias factor in this simple model, 
Eq. (7), can be shown to change by a factor of about two 
when the the mean flux, Eq. (5), changes within its 2-a 
interval, similar to the variations seen in Fig. 17 of Croft 
et al. (2000). Equivalently figure 10c of McDonald (2001) 
shows that a 3% change in the mean flux can lead to a 
40% change in the power spectrum amplitude. 

As pointed out above, the validity of Eq. (7) is not re- 
stricted to the linear regime, since there are other small 
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parameters that enter. The value of Aa^ is eontroUed by 
the mean flux as given by Eq. (5). Since the mean flux 
has to be between and 1, in particular F 0.7, Aa^ is 
constrained to be smaller than unity^, even in the limit 
that cr^ ^ 1. Eq. (7) then follows as long as ^ <C cr^, irre- 
spective of the value of ^ (since Aa^ 1, <C 1 in this 

regime). The condition that ^ ^ cr^ means in the Ly-a 

case that the scales under consideration be much larger 
than the thermal broadening scale, k <C fcxH, the charac- 
teristic smoothing scale where is evahiated. This is in 
any case where most of the information can be extracted, 
so the approximations leading to Eq. (7) are physically 
relevant. 

Of course, when ^ approaches unity, the Gaussian model 
for the density field PDF is not a good approximation any- 
more. In order to see how the arguments above generalize 
to the non-Gaussian case, let's consider the expansion of 
the two-point PDF when ^ <C <t^, 

P(^i, ^2) ~ P{Si)P{S2) [1 + biSi)bi52) (8) 

where b{6) = S/a"^ in the Gaussian case. Although the 
relation in Eq. (8) was written starting from a Gaussian 
PDF, its validity is much more general. It can be ob- 
tained rigorously in the weakly non-linear regime (in tree- 
level PT, Bernardeau 1996), and in the highly non-linear 
regime for factorizable hierarchical models (Bernardeau & 
Schaeffer 1999), where specific predictions for the function 
b{S) can be given in both cases. Then, the linear relation 
in Eq. (7) between S^p and ^ continues to hold, but with a 
different bias factor than in the Gaussian case. Note that 
in the hierarchical model, ^ is understood to be the fully 
non-linear correlation function. 

These two examples, however, are not exactly relevant 
for structure formation at z ^ 3 at the scales of inter- 
est. Tree-level PT is not a good approximation because 
the spectral index at the non-linear scale is so negative 
i^cS ^ —2.5 that loop corrections to moments of the den- 
sity field are very important, as noted above for the power 
spectrum. On the other hand, we cannot assume a factor- 
izable hierarchical model, which is expected to become a 
good approximation only in the highly non-linear regime. 

However, the structure of loop corrections still respects 
Eq. (8) in the following broad sense. Consider the calcula- 
tion of connected two-point moments, ( (5f(5| )c, which fully 
characterize the connected two-point PDF. At large scales 
it follows that (Bernardeau 1996) 

(<5f<5|),^Cp,af+«-^)a, (9) 

where Cpq are numbers that depend weakly on scale 
through derivatives of the linear variance and smoothed 
two-point function with respect to scale {L denotes linear 
quantities). Since all connected moments are proportional 
to the connected two-point PDF should also be pro- 
portional to ^L, as in Eq. (8), where ^ is replaced by ^l- In 
addition, it can be shown that Cpq = CpiCqi, which leads 
directly to the factorization b{5i)b{52), with b{S) calcula- 
ble from tree-level PT. As loop corrections are included, a 
subset of loop diagrams will lead to non- linear corrections 

^This can breaJj down for non-Gaussian distributions, since as 
a increases the PDF becomes more peaked about 5 = — 1 and one 
can thus tolerate a larger A. For the models discussed in Section 2, 
Aa^ » 1. 



to and ct^ in Eq. (9), making £,l ^ £, and (t| in 
Eq. (9). In addition, there will be new links in the dia- 
grams leading to higher powers of and u\\ however, at a 
given order in PT since ^i, a\ the latter are much more 
important than the former. As a result of this, we expect 
that as non- linear scales are probed, as long as ^ <C cr^, 
Eq. (9) becomes 

((5[5|),~Cp,(a)a2(^'+''-2) (10) 

where we have absorbed additional a dependence into 
Cpq. Note that this also respects the scaling expected in 
the highly non- linear regime, as long as Cpq{a) const, 
as (J ^ oo. Recent studies of the behavior of (^S^i&\)c as 
a function of a and ^ in numerical simulation support the 
above arguments for ^ <C cr^ (Gaztanaga, Fosalba & Croft 
2001). However, one must keep in mind that they studied 
CDM models at z = rather than z ~ 3 as concerns us 
here'''. 

Now let's consider the effect of redshift distortions. In 
this case, the redshift-space position s = x— /^/^(x), where 
/ « « 1 at the redshifts of interest and is the line of 
sight velocity, normalized by the Hubble constant so that 
in linear PT, V • u = 5. The transformation between opti- 
cal depth in real and redshift space reads Tgds = rdz (e.g. 
Hui et al. 1997b ; Gaztafiaga & Croft 1999) with ds = Jdz 
and ,/ = |1 — fVzUzl- This leads to a three- dimensional 
fiux power spectrum (following the calculation in Scocci- 
marro, Couchman & Prieman 1999), 

Pf(k) = yrfVe*'-(e'^--f^"-JiJ2e-^i/-^ie-^=/^^), (11) 

where r = xi — X2 and Auz = Uz{xi) — Uz{x.2) . At large 
scales, one is allowed to expand all the exponential factors 

and use 5,V 1, in which case e'^^'^^^' J1J2 » 1 + 
0{\/zUzf, and the remaining exponentials give Pf{k) = 
{Ai^)'^{l + /Sji'^fPLik) (Hui 1999), where we assumed r = 
A{1 + 6y and = kz/k, [3 = f /u. This is analogous 
to the standard result for the galaxy redshift-space power 
spectrum in the linear regime (Kaiser 1987). 

Non-linear corrections break the relation 6 = V ■ u = 6, 
valid in linear PT. In particular, a similar calculation to 
that in Scoccimarro & Frieman (1996) for the density field 
shows that the critical index for the velocity divergence 
power spectrum is « —2. Therefore, for the models un- 
der consideration with n « — 2 one has ^ = V- u~Jz, < 5; 
velocity fields are much more "linear" than the density 
field. We see from Eq. (11) that at small scales the J fac- 
tors should suppress the regions undergoing turnaround, 
where J ~ 0. These regions, however, correspond to 
1 = VzMz ~ 0/2, ~ 5z,/3, thus 5l ~ 3, which corresponds 
to very high values of 5 and were suppressed already in 
real space by the exponential dependence on density. On 
the other hand, the effects of the pairwise velocity along 
the line of sight encoded in g**^^^^"^ plays a role, at high 
enough kz oscillations of this factor damp power along the 
line of sight leading to a negative quadrupole to monopole 

*In any case, even if there were higher-order corrections in ^ in 
Eq. (10). still what enters is the non-linear ^, as shown in Scoc- 
ciiiiarro (2001). Therefore, our arguments about the difficulties of 
determining the initial shape of the dark matter power spectrum 
from the forest still apply. 
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ratio, as in the case of the density field. However, since 
velocity fields should also display critical behavior as den- 
sity fields, redshift space corrections should become very 
similar for models with different initial spectra. Therefore, 
we believe that the effects of velocities do not affect the 
conclusions above. 

5. CONCLUSION 

We have shown that at 2; ~ 3 the constraints on the pri- 
mordial power spectrum of density perturbations that can 
be obtained from the Lyman alpha forest flux power spec- 
trum are significantly affected by non-linear corrections to 
the evolution of the density perturbations. Because the 
primordial index of the power spectrum is so negative on 
the scales of interest non- linear corrections are much larger 
than naively expected. These corrections drive the power 
spectrum towards a critical slope, n ~ — 1.4 irrespective of 
the primordial spectral index. This makes the determina- 
tion of the primordial spectral index from measurements 
of the flux power spectrum particularly difficult. 

If the objective is to constraint the slope of the pri- 
mordial power spectrum, over a wide range of models in 
the relevant part of parameter space non-linear corrections 
are very large. Figure 4 shows that all the scales probed 
by the forest are significantly affected by the non-linear 
corrections for some shapes of the power spectrum. Thus, 
when trying to constrain cosmological parameters, the flux 
power spectrum cannot be used as a probe of the linear 
power spectrum of fluctuations at z ~ 3 except at the 
largest scales. Figure 1 illustrates how this "degeneracy" 
manifests itself when determining rip, and figure 2 shows 
an example of two models with originally very different 
Tip that become almost indistinguishable due to non-linear 
evolution. This problem is made worse by the sensitivity 
of the bias to F. 

There are several ways to increase the sensitivity of the 
forest to the shape of the primordial power spectrum. First 
future improvements in measurements of the equation of 
state might help reduce the error on the shape. Second, 
the Sloan Digital Sky Survey (SDSS) has the capability 
of extending flux power spectrum measurements out to 
larger (and therefore more linear) scales, which will help 
towards tightening constraints on the shape of the power 
spectrum. 
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